Purification through Zeno-like Measurements 
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A series of frequent measurements on a quantum system (Zeno-like measurements) is shown to 
result in the "purification" of another quantum system in interaction with the former. Even though 
the measurements are performed on the former system, their effect drives the latter into a pure 
state, irrespectively of its initial (mixed) state, provided certain conditions are satisfied. 

PACS numbers: 03.65.Xp 
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It is well known that unstable particles or quantum 
states display a peculiar behavior at short and long 
times (JJ. Phenomenologically, they are known to de- 
cay exponentially and this is well confirmed experimen- 
tally [2j. Short-time deviations were only observed very 
recently Q. The deviations from the familiar exponen- 
tial decay law are unavoidable consequences of the quan- 
tum dynamics both at short and long times, and the 
derivation of the exponential decay law itself is not a 
trivial matter in quantum mechanics. These deviations 
reflect the unitarity of the time evolution operator or the 
time reversal symmetry of the Schrddinger equation at 
short times and the lower boundedness of the Hamilto- 
nian or the stability of the vacuum at long times. See, 
e.g., Ref. [H for a review. 

The quantum behavior of unstable states at short times 
has been one of the central issues of investigation and dis- 
cussion in recent years, since it is closely connected to the 
so-called quantum Zeno effect (QZE) 0,0, where the act 
of measurement @ (usually represent^ the von Neu- 
mann projection or the generalized spectral decomposi- 
tion |2j) affects in an essential way the dynamics of the 
measured system and results in a hindrance of the decay 
process. The first attempt at the experimental observa- 
tion of the QZE in an atomic transition process ||, fol- 
lowing Cook's theoretical work has triggered heated 
discussions on this subject. Furthermore, another excit- 
ing experiment has been reported very recently: the ob- 
servation of the QZE (and also of the inverse QZE [l(ij]) 
in an atomic tunneling process [ill ] , which is the first ex- 
perimental observation of the (inverse) QZE in a truly 
unstable quantum system, unlike in the previous experi- 
ment 8] performed on an oscillating system. 

In this Letter, we will shed new light on another (and 
so far not well explored) feature of the quantum dynam- 
ics with measurements, closely related to the QZE. No- 
tice first that the system under consideration cannot be 
considered completely isolated and usually interacts with 
other systems. Therefore, it would be interesting and 
maybe more realistic to consider the case where the mea- 
surement, represented by a von Neumann projection for 
simplicity, is not performed on the total system, but only 
on the system of interest. Here, we consider such mea- 
surements and address the following point: How does a 



series of frequent measurements on a system affect the 
dynamics of another system in interaction with the for- 
mer? Under frequent measurements performed only on 
the former system, the latter evolves away from its initial 
state. We shall show that such measurements can result 
in a "purification" phenomenon. That is, a series of fre- 
quent measurements on system A, represented by pro- 
jections on a given (usually pure) state of A, makes the 
state of system B, which interacts with A and is initially 
in any (mixed) state, approach a pure state, if certain 
conditions prescribed below are satisfied. 

Let a total quantum system A+B be described by a 
Hamiltonian H of the form 



H = H A + H E 



(1) 



where H^b) stands for a free Hamiltonian of system 
A(B) and H mt for an interaction. We prepare system 
A in its initial (pure) state |0)(0| at time t = 0. The ini- 
tial state of system B, denoted by pb, can be arbitrary. 
The initial state of the total system is 



(2) 



and its dynamics is governed by the Hamiltonian Q un- 
less it is interrupted by a series of measurements on sys- 
tem A, each of which is represented by a projection oper- 
ator O = |0) (0| ® 1, performed at time intervals r. No- 
tice that this operator, even if it is a bona fide projection 
operator (we assume that |0) is normalized), does not 
return the time-evolved total system to its initial state. 
The projection is partial, in the sense that only the state 
of system A is set back to its initial state and that of 
system B is not initialized, even though the dynamics of 
B is certainly affected. 

After TV such measurements have been done, the sur- 
vival probability of finding system A still in its initial 
state is represented by 



P (T \N) 



Tr 
Tr B 



(Oe- lHr O) N p (Oe iHr O) N 

iHr I i 



(3) 



Notice that the quantity V^(r) = (cp\e '"' \<p) is an oper- 
ator acting on the Hilbert space of system B. The density 
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operators of the total and B systems read 

p iT) (N) = {Oe- lHT O) N p (Oe lHT O) N /P (T) (N) 

= |0>^| ®/4 T) (AO, (4a) 

P [ b\n) = (V,(t)) N Pb (vI(t)) N /PV(N), (4b) 

respectively. We only collect the right outcomes of mea- 
surements: this is implicit in the normalization factors 
in |4"|. Experimentally, this means that after each mea- 
surement, only those events will be retained in which 
system A has been found in its initial state. 

In the ordinary situation, one performs infinitely fre- 
quent measurements by taking N — > oo and r — > 0, 
keeping Nt = T, a finite nontrivial value; one easily 
checks that the ordinary QZE [l| appears in this case 
and the survival probability p( T \N) increases as N be- 
comes large, approaching unity in the N — > oo limit [5|- 
At the same time, the dynamics of system B becomes uni- 
tary in this limit, and this is an example of the so-called 
"quantum Zeno dynamics" jl2l | . However, we stress that 
our interest lies in a different situation: we keep the time 
interval r between measurements finite and nonvanish- 
ing. If N were taken to be oo, the survival probability 
p( T )(jV) would decay out completely for such but 
we are interested in the asymptotic behavior of the state 
of system B for large but finite values of N. We ex- 
pect that the effect of repeated measurements on system 
A would modify the dynamics of system B through its 
interaction with the measured system A, even if B has 
never been directly measured. To examine this idea, we 
need to clarify the asymptotic behavior of the state of 
system B, pg(N), for large N. 

It is clear that the behavior of p%\N) is governed by 
the operator V^(r) in (@J. Let us consider its eigenvalue 
problem. Since this operator is not hermitian, VJ(t) ^ 
V^r), in general, we need to set up both the right- and 
left-eigenvalue problems 

V^(t)|u„) = \ n \Un), (v n \V<f,(T) = X„(v n \. (5) 

The eigenvalue A n is in general complex-valued. Let us 
assume that the spectrum of the operator V$(t) is dis- 
crete and nondegenerate, and its eigenvectors form an 
orthonormal complete set in the following sense 

^ \u n )(v n \ = 1, (v n \u m ) = 5 nm - (6) 

n 

[It will soon become clear that the assumption of the 
nondegenerate spectrum is not essential for the follow- 
ing discussion except for that of the largest (in magni- 
tude) eigenvalue A max .] The operator itself is expanded 
in terms of its eigenvectors 

V^T) =J2^nK)(v n \ (7) 



and we obtain 

n 

One can show that the absolute value of the eigen- 
value A n satisfies the inequality < |A„| < 1, Vn, which 
reflects the unitarity of the time-evolution operator. It is 
now evident that, in the large N limit, the operator (JHJ 
is dominated by a single term 

(V<p(t)) N ' argC N > A^ ax |ti max )(u ma x| ! (9) 

where |u max ) and (i> max | are the eigenvectors belonging 
to A max , provided the largest (in magnitude) eigenvalue 
A m ax is unique, discrete, and nondegenerate. 

Thus we reach the conclusion, under the assumption 
of unique, discrete, and nondegenerate A ma x, that, in the 
large N limit with a nonvanishing r, the state of system 
B in interaction with system A, on which N measure- 
ments are performed at time intervals r, asymptotically 
approaches the pure state |iimax) 

(t)/tv7-\ large N , , . . , , ,. A , 

P B [N) > |Uma X )(Wma X |/(Wmax|u ma x), (10) 

with probability 

p( T )(j\T) largC N > |A m a X | 2Ar (Wma X |Umax)(Wma X |PB|Wmax)- 

(11) 

Notice that the final pure state |u max ) is independent 
of the choice of the initial state of system B, i.e., any 
initial (mixed) state shall be driven to the unique pure 
state |iimax) by repeated measurements performed on the 
other system A. Since the asymptotic state |u m ax) is one 
of the eigenstates of the operator (r) , we have the pos- 
sibility of adjusting the interaction strength and the mea- 
surement interval and of choosing an appropriate initial 
state \(j>) for system A so that a desired pure state |u m a X ) 
is realized in system B after a large number of measure- 
ments on A [as long as the probability (N) does not 
become meaninglessly small] . This discloses another fea- 
ture of the Zeno phenomenon: the action of the quantum 
Zeno-like measurements dramatically affects the dynam- 
ics of system B of interest. 

A few comments are in order. First, the existence of 
a unique, discrete, and nondegenerate A max of the op- 
erator V^,(t), which has been assumed here, is essential 
for this purification mechanism. Even though this con- 
dition is satisfied for some systems with a discrete spec- 
trum as will be shown in the examples below, some defi- 
nite mathematical criteria for its validity have yet to be 
clarified. In particular, it remains open if the present 
purification mechanism can be applied to systems with 
continuous spectra. Nevertheless at the same time, it 
would be worth stressing that not a few discrete quan- 
tum systems, including 2- or 3-level systems which play 
important roles in the field of quantum information and 
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computation, certainly fall into the category of systems 
with unique, discrete, and nondegenerate A max . Second, 
it is easy to see that the approach to the final pure state 
l^max) is governed by the ratio between the largest and 
the second largest (in magnitude) eigenvalues of the op- 
erator (t) . It is possible that as the number of degrees 
of freedom increases, the eigenvalues A„ distribute more 
closely to each other, which would make the present pu- 
rification process less effective. Lastly, even though the 
measurements are repeated many times, like in the bona 
fide Zeno case, the present scheme does not explicitly 
rely on the peculiar quadratic behavior of quantum sys- 
tems at short times. As far as the essential assumption 
on the spectrum of the operator V^(r) is satisfied, there 
will be no limit on the time interval t. Notice that the 
repetition of one and the same quantum measurement is 
crucial here. 

Let us illustrate the above conclusion in a simple but 
still nontrivial model. We consider two single-mode har- 
monic oscillators a and b, in interaction in the rotating- 
wave approximation. The total Hamiltonian reads 



H = Via) a + vtfb + ig{a)b - afe 1 *), 



(12) 



where the frequencies VL and u> and the coupling constant 
g are real parameters. The spectrum is discrete. We 
prepare system A (oscillator a) in some definite pure state 
\(j>) (typically a number state \n a ) or a coherent state 
a)) at time t = and let it evolve under the above 
Hamiltonian. Then the initial state of system A starts 
to evolve towards other states owing to the coupling to 
system B (oscillator b), the initial state of which can be 
arbitrary. The state of oscillator a is projected onto its 
initial state \<p) at each measurement, and the interval 
between measurements r is taken small, compared with 
the typical time scales of the system, e.g., 2tt/S in 114|) 
below. 

The eigenvalue problem (JSJ of the relevant operator 
V^(r) is solved exactly in this case. Indeed, since the 
time-evolution operator e~ lIiT can be factorized as 



iHi 



e Aa t b e Ba f a e Cb 1i b e -Aab'' 



in terms of the r-dependent functions 
{g/5) sin St 



A = 



cosSt + i[(Q — uj)/2S] sin<Sr ' 



B = --(Q, + w)t - In 



C= --(n + w)r + ln 



. .fi — U) . r 
cos or + i — — — sin or 
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cos St + i - 



to 



2S 



sin St 



(13) 



(14a) 



, (14b) 



(14c) 



where S = \J g 2 + (£1 — w) 2 /4, we easily find the eigenvec- 
tors \u n ) and (v n \ of the operator V^(r), once the initial 
state |0) of oscillator a is specified. 



If we prepare oscillator a in the number state \n a ) at 
t = 0, the relevant operator is calculated to be 



I' a 



Vnjr) 



II, 



2 [(n a -fc)!] 2 fc!' 



IB 



{ _ A 2 e C )na 



x e 



n a —k 

' n ^ b+e ^ 



(15) 



from which we understand that the number states jn?,) 
of oscillator b constitute the set of eigenvectors of the 
operator l|15l) . Therefore the state of oscillator b is driven 
to a number state irrespectively of its initial state, when 
the coupled oscillator a is repeatedly confirmed to be in 
the number state \n a ). The state of system B is purified 
into a number state. 

On the other hand, when oscillator a is prepared in a 
coherent state \a) and found to be in this state at every 
r, the relevant operator is rearranged to be 



)]M 2 e D(&+,6) 



V a (T)=e-^- e - A2 ^- 
where the operator D(b\ b) is expressed as 



D(b\b) = C 



Aa* 



1 



b- 



Aa 



1-e- 



(16) 



(17) 



It is easily understood that the state of oscillator b ap- 
proaches a coherent state b\f3) = j3\(3) with (3 = Aa/(1 — 
e~ c ), since this is the right-eigenvector of D(b^,b) be- 
longing to zero eigenvalue and therefore that of V„ (r) 
belonging to the largest (in magnitude) eigenvalue |14| . 
The state of system B is again purified into (another) 
pure state \(3). 

In Fig.^ the survival probability P^(N) and the so- 
called fidelity 



F^(N) = (wIp^CaOK 



)/(« 

max | ^max ) (18) 



are shown as functions of the number of measurements 
N for the case ljl7>)l - lfr7|) with a particular choice of pa- 
rameters. In order to make the purification procedure 
more effective, it is preferable that (i) the magnitude of 
the largest eigenvalue of the operator V^(r) be close to 
one, |A max | ~ 1, which maintains the probability P( T >(N) 
large enough even for large TV [see (jllf) ] , and (ii) the other 
eigenvalues be all small (in magnitude) compared with 
I A m ax |, in order to realize a faster approach to the final 
pure state |u m ax). For this purpose, one may adjust the 
relevant parameters, such as the interval between mea- 
surements r, the strength of the interaction g, and the 
state | a) onto which system A is projected. The condi- 
tion | A max | — 1 is satisfied in general if the interval t is 
taken to be small enough as in the ordinary Zeno mea- 
surements, but one can optimize this procedure and find 
better values of r, not necessarily very small .13], that 
satisfy both conditions (i) and (ii). See Fig.QJand its cap- 
tion, where r is tuned so as to satisfy the conditions (i) 
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FIG. 1: Probability P {t) (N) and fidelity F (r) (iV) for the 
model 11 1 2> when the initial state of system A, onto which N 
measurements are performed, is a coherent state \a) and that 



of system B is ps oc e 



ub^b/k B T 



(thermal, i.e., maximally 



mixed). The parameters are taken to be u) = 1, g = 0.2, 
T = 1, a = 0.5 and r = 27r/[(fi + w)/2 + 5] ~ 5.24 in units 
such that 7i = fes = f2 = 1. ris tuned so as to satisfy 
the condition |A max | = 1 13], and the ratio of the second 
largest (in magnitude) eigenvalue to the largest one A max is 
\e c \ = yi-( 5 /5) 2 sin 2 5r~0.37 0. 



and (ii) for the case (|16fl - l|17[) . and the purification mech- 
anism becomes very effective after only N = 2 steps. 

The above arguments clearly and explicitly show how 
the action of repeated measurements (projections) on one 
system A can affect the dynamics of the other system B 
in interaction with A. Interestingly enough, even though 
the effect of the measurement on the latter system B is in- 
direct, its influence is far-reaching if the measurement is 
repeated many times: irrespectively of its initial (mixed) 
state, the state of system B is purified towards a pure 
state, provided the conditions on the spectrum of V^(r) 
are satisfied. The final state of system B is prescribed by 
the total Hamiltonian, the pure state (usually taken to be 
the initial state) onto which system A is projected by the 
measurement, and the time interval between successive 
measurements. This opens a new possibility on how to 
control the state of a quantum system on which we have 
no direct access. If another system under our control can 
be coupled to the former system, we would only have to 
decide which state has to be measured on the control- 
lable system. After such measurements are performed 
many times at the prescribed time intervals, the desired 
pure state would be realized with some probability in the 
system beyond our control. 

Purification of quantum states is now considered to 
be one of the key technologies for quantum information 
and computation, and is being widely explored |l5j (es- 
pecially in the context of "entanglement purification"). 
Compared to some other procedures, the idea here is 
rather simple: one has only to repeat the same mea- 



surements. The objects to which the present method 
is applicable are general and not restricted to "qubits." 
Furthermore, it is worth emphasizing the versatility of 
this procedure, i.e., the possibility to adjust the target 
pure state, the balance between fidelity and probabil- 
ity yield, and so on. These issues would deserve further 
study, for example, in the context of quantum informa- 
tion and computation. 
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